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ABSTRACT 

In this study, a laboratory soil column experiment was 
carried out for estimating the dispersivity (A,) in sandy 
loam soil by the inversion of the Advection- 
Dispersion Equation. An 89 cm deep and 15 cm 
diameter cylindrical column was filled with sandy 
loam soil in such a way that the bulk density of the 
soil in field condition was strictly satisfied in the 
column. Potassium bromide was used as the 
conservative tracer in the study. The tracer 
concentration was allowed to pass through the column 
at a constant rate and the leachate collected at 
different time periods were tested for the 
concentration of bromide. The experiment was carried 
out until the relative concentration (c/cO) has reached 
1. The calculated relative concentration was plotted 
against time for obtaining the observed breakthrough 
curve (BTC). Another BTC was plotted using the 
analytical solution of Advection Dispersion Equation. 
The effective dispersion coefficient (De)value which 
gave the best fit between these two curves was 
chosen by trial and error method and the dispersivity 
of the soil was calculated (De = Av)for a constant pore 
water velocity (v). The dispersivity of the sandy loam 
soil under study was found to be 1.7 cm which 
matches well with the previous researches. 

Keywords: Dispersivity: Advection Dispersion 

Equation, Effective dispersion coefficient, Laboratory 
soil column. 

INTRODUCTION 

Understanding solute transport processes between the 
soil surface and the groundwater table is essential for 
limiting or effectively managing soil and groundwater 
pollution. Dispersivity are typically derived either 
from observed depth profiles of inert tracer 


concentrations or from breakthrough curves (BTC) of 
tracer concentrations that are measured in the effluent 
of columns or that are measured in a soil profile using 
solution samplers or other devices such as time 
domain reflectometry (TDR) probes that monitor 
solute concentrations. 

The dispersivity (X) reflects the degree of mechanical 
mixing of solutes, caused by variations in the pore 
water velocity. The soil-column experiment was 
widely used to evaluate the transport model and 
determine the fate and migration of contaminants 
through soils. The value of X has traditionally been 
considered under saturated water condition and 
usually reported in the range of 0.1-2 cm for 
homogeneous saturated soils (Bear, 1972).Silliman et 
al., (1987) studied the scale effect in dispersion of 
solutes in porous media under laboratory conditions 
and the results for various sand-packing arrangements 
revealed that the breakthrough curves for uniform 
coarse sand showed a constant dispersivity of 
approximately 0.02 m. Huang, K. et al., (1995) 
conducted an experimental investigation of solute 
transport in large, homogeneous and heterogeneous, 
saturated soil columns. Values of the dispersivity (X) 
for the homogeneous sandy column ranged from 0.1 
to 5.0 cm, while those for the heterogeneous column 
were as high as 200 cm. 

MATERIALS AND METHODS 

Laboratory Soil Column Experiment 

Tracer (Bromide) movement through the saturated 
agricultural soil was investigated in soil columns. The 
dispersivity of bromide ions in sandy loam soil was 
measured in this study. Potassium bromide was used 
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as the conservative tracer. Bromide salts have proven 
to serve as excellent tracer materials since they are 
highly soluble, non-degradable, relatively 
inexpensive, and can be measured quantitatively in 
very low concentrations. 

A cylindrical laboratory soil column of 89 cm height 
and 15 cm diameter (Fig. 1) was used to find the 
dispersion of potassium bromide tracer in a 
homogeneously packed sandy loam soil. The bulk 
density of the soil measured using core cutter method 
was found to be 1.58 g cm' 3 . The soil column was 
filled in such a way that the bulk density of the soil in 
the field was strictly satisfied in the laboratory soil 
column. 


Porosity = 


Bulk Density 
Particle Density 


( 4 ) 


The uniform soil column was saturated initially by 
passing water through it. The potassium bromide 
(KBr) solution having an initial concentration of 20 g 
per 200 ml of tap water was allowed to pass through 
the soil column from its upper end and the leachate 
was collected at the lower end. The input solution was 
given at a constant rate of 800 cm h' throughout the 
experiment. Since a constant input of solute at a given 
concentration (cq) was given at the surface boundary, 
this is known as step input. Concentration of bromide 
ions in the leachate collected at different time steps 
were determined by Mohr’s (Argentometric titration) 
method. 



Volume of the soil column, 

V = nr 2 h ( 1 ) 

= ttx(7.5) 2 x(89) = 15719.63 cm 3 
Bulk Density of the soil ( BD ) = ^ (2) 

where, M is the mass of the soil, g and, V is the 
volume of soil, cm 3 . 


Mass of the soil to be filled in the soil column, 
M = BD x V (3) 

=1.58 x 15719.63 
= 24837.01 g =24.84 kg 
i.e. for satisfying the bulk density of the soil (1.58 g 
cm' 3 ), 24.84 kg of soil was filled in the soil column 
providing an additional free board of 1 cm at the top. 
The porosity of the soil was found to be 0.37 (Eq. 4). 


Mechanism of solute transport in soils 

'The solute conservation equation states that the 
change in total solute concentration (Cf) in units of 
mass per volume of soil in a representative elementary 
volume over time is equal to the negative of the rate 
of change of the solute flux with depth, minus any 
source/ sink ( S) in units of mass per volume per time 
(Radcliff and Simunek, 2010). Only solute flow in the 
vertical direction (z positive up) is considered for 
simplification purposes. 

(s\ 

dt dz ^ ' 


where, Cjis the total solute concentration in units of 
mass per volume of soil, S is the source/ sinks in units 
of mass per volume per time and,/ s is the total solute 
flux. 


This is the partial differential equation (PDE) with 
independent variables z and t and the dependent 
variable C T . The units of each term are [ML 3 T 1 ]. 
Solute flux (Js) i s the rate at which the solute flows 
through a given cross-sectional area. 

The total solute concentration (Cf) consists of the 
dissolved or liquid-phase solute concentration (c), 
expressed as a mass of solute per volume of liquid, 
and the adsorbed or solid-phase solute concentration 
(s), expresses as a mass of solute per mass of soil. To 
convert each of these concentrations to units of mass 
per volume of soil, c is multiplied by 6 and s is 
multiplied by p b . 

C T = 9c + p b s ( 6 ) 

where, c is theliquid-phase solute concentration [ML' 
3 ], s is the solid-phase solute concentration [MM 1 ], 0 
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is the volumetric water content [L 3 L“ 3 ] and ,p b is the where, A is the dispersivity, [L] and, v is the pore 
bulk density of soil [ML' 3 ]. water velocity, [LT 1 ]. 


Solute flux in the liquid consists of three components 
that represent different transport processes: advective 
(also called convective) flux (J Ic ), diffusive flux (Ju) 
and hydrodynamic dispersive flux (JjA. 

Is = Jic + Jid + Jih ( 7 ) 

Advection Dispersion (Convection- Dispersion) 
Equation 


Most mechanistic transport models for solutes in 
porous materials are based on the convection- 
dispersion equation (CDE). The CDE is a partial 
differential equation representing mass continuity for 
movement of a given solute in a porous medium by 
dispersion and convection under specified initial and 
boundary conditions. Appropriate terms are 
incorporated into to the CDE to account for chemical 
and/or physical sink/source interactions between the 
solute and porous medium. The CDE can be 
developed microscopically based on Brownian 
motion, or macroscopically, based on Fick’s law. 

The final solution to the Advection-Dispersion 
equation obtained by the Laplace transform method 
is: 


c(z,t) 


= c 0 


v 2 t 


nD, 


•exp | 


(z — vty 
4 ~D~t 


+ 2 er f c I 


z — vt 

vy/4 ~Dgt/ 


1 vz 

-2\ 1 + K 


V 2 t\ (VZ\ 

+ d7 exp U er/c 


z + vt 
4D e t; 


( 8 ) 


where, erfc is the complementary error function 
defined as er/c(x) =-p/^e - * 2 dx, c is the 

concentration at depth z at time t [ML 3 ], Co is the 
initial concentration [ML' 3 ], D e is the effective 
dispersion coefficient [L T" ], v is the pore water 
velocity [LT 1 ], z is the depth [L] and, t is the time 
[T], 

Effective dispersion coefficient, 

D e = Av 

( 9 ) 


Determination of Bromide Ion Concentration by 
Argentometric Titration (Mohr’s method) 

Mohr’s method determines the bromide ion 
concentration of a solution by titration with silver 
nitrate (AgN0 3 ). As the silver nitrate solution is 
slowly added, a precipitate of silver bromide forms. 

Ag (aq) T Br ( a q) » AgBr( s ) 

The end point of the titration occurs when all the 
bromide ions were precipitated. Then additional silver 
ions react with the chromate ions of the indicator, 
potassium chromate, to form a red-brown precipitate 

of silver chromate. 

+ 2 - 

2 Ag (aq)"*" Cr04 ( a q) > Ag 2 Cr 04 ( s ) 

The pH of the sample solutions should be between 7 
and 10 because chromate ion is the conjugate base of 
the weak chromic acid. A suitable pH was achieved 
by saturating the analyte solution with sodium 
hydrogen carbonate (NaHC0 3 ). 

Although the silver bromide that forms is a white 
precipitate, the chromate indicator initially gives the 
cloudy solution a faint lemon-yellow colour. The 
endpoint of the titration was identified as the first 
appearance of a red-brown colour of silver chromate. 
The titration was repeated with further aliquots of 
diluted samples until concordant results were 
obtained. 

A blank titration is carried out by titrating a fixed and 
known concentration of titrant into a solvent with zero 
analyte. The only difference from the regular titration 
is the absence of analyte. This allows the amount of 
reactive substance within the plain solvent to be 
determined and hence allows a determination of the 
error in future titration experiments using this 
solvent.The blank titration was conducted to correct 
the errors that might occur in the titre values due to 
the presence of chloride impurities in the potassium 
chromate indicator.The volume of AgN0 3 needed for 
the complete precipitation of the bromide ions (titre 
value) was noted and the concentration of bromide 
ions in the leachate was calculated using the following 
equation: 

(Vs~Vb) x M AgN03 x Equivalent wt of Br xlOOO 
^ Vol. of sample taken (^9) 
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where, c is the concentration of bromide in the 
leachate in mg l' 1 , V s is the volume of AgNC >3 needed 
for the sample in ml, V B is the volume of AgNC >3 
needed for the blank in ml (0.2 ml) and,M Ag N03 is 
themolarity of AgNC >3 (0.1M). The equivalent weight 
of bromide is 79.9 g. 

The concentration (c) of the leachate and the relative 
concentrations ( c/co ) for different time steps were 
calculated using Eq. 10. Once the concentration of the 
leachate was calculated, a graph was plotted with the 
relative ( c/co ) against the time (t). This graph is said to 
be the observed breakthrough curve. 

Break Through Curves (BTCs) 

The concentrations can be plotted as a function of 
time at a given depth for various values of effective 
dispersion coefficient ( D e ) and the plots are called 
breakthrough curves (BTC). Break through curves are 
extensively used to characterize the physicochemical 
processes involved in the transport of solutes in 
porous media. Analysis of column BTCs can provide 
useful estimates of the parameters for the 
physicochemical processes involved in subsurface 
solute transport (Radcliff and Simunek, 2010). Such 
parameters are essential for developing and validating 
theoretical models for solute transport. 

The relative concentration (c/co) obtained from the 
Advection-Dispersion Equation (Eq. 8) was plotted 
against corresponding time periods (t) to generate the 
BTC (predicted) at 89 cm depth. Microsoft excel was 
effectively utilised to calculate the relative 
concentration values. Excel allows the user to 
calculate the complementary error function directly. 
The effective dispersion coefficient is a constant, 
rather than a function of time or mean travel distance 
(Dagan, 1988; Neuman and Zhang, 1990; Rajaram 
and Gelhar, 1993). The effective dispersion 
coefficient (D e ) in Eq. (8) was chosen arbitrarily by 
trial and error method. 

The pore water velocity (v) was calculated as: 


where, 0 is the porosity and,/^ is the Darcy velocity. 

Jw = 7 ( 12 ) 


the soil column in cm 2 . The pore water velocity(v) 
was calculated as 12.24 cm h' 1 (i.e. 293.79 cm d' 1 ). 

The BTC for observed bromide concentration was 
obtained using Eq. 10 and the resulting plot was fitted 
to the solution of 1-D ADE (Eq. 8). The effective 
dispersion coefficient (. D e ) which gives the best fit 
between the observed and predicted BTC’s was 
chosen by trial and error method. Finally, the 
dispersivity (2) for the bromide ions in the sandy loam 
soil column at a depth of 89 cm was calculated using 
Eq. 9. 

RESULTS AND DISCUSSIONS 

The hydrodynamic dispersion coefficient (. D e ) and the 
pore water velocity (v) were derived from depth 
profile (89 cm) assuming that v and D e are constant in 
the soil profile and do not change with depth and 
time, i.e. a hydro dynamically homogeneous soil 
profile and steady state flow conditions. The 
analytical solutions of the solute transport equation 
(ADE) was derived and fitted to the observed BTCs. 



Fig. 2: Relative Concentration vs. Time curve 
obtained from the estimated bromide 
concentration 


An initial concentration of 20 g per 200 ml of 
potassium bromide was given as the input tracer 
solution. Initially, the concentration of bromide was 0 
in the leachate and as time passes the bromide starts to 
appear in the leachate. The first appearance of 
bromide was found 43 minutes after the starting of the 
experiment. The concentration of bromide in the 
leachate increased with time and reached 100 % 
concentration after 12.25 hrs. 


where, Q is the constant discharge rate of the step The bromide concentration of the leachate (c) 
input (800 cm 3 h _1 ) and A is the input surface area of estimated using Mohr’s test was determined by Eq. 

(10) and the relative concentration (c/cq) for different 
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time steps were calculated. The observed BTC (Fig. 
2) in the homogeneous sandy loam column was found 
to be a smooth and sigmoidal (S-shaped) curve. After 
about 6.67 hrs, the relative concentration reached a 
value of 0.5 and half of the BTC had arrived. 
Eventually, concentration reached the input 
concentration at 12.25 hrs. At this point, solute has 
arrived at this depth even in the smallest pores. 

The effective dispersion coefficient ( D e ) and 
dispersivity (2) were determined by fitting the 
analytical solution of the convection-dispersion 
equation (Eq. 8) to the BTC obtained by estimated 
relative concentration. The solution of 1-D ADE was 
fitted to the locally measured BTC by randomly 
changing the value of effective dispersion coefficient 
(. D e ) by trial and error method. Since the dispersion of 
the solute front was affected by D e , so was the BTC. 
As D e increased, the BTC became more dispersed. 


♦ c/c 0 (e sti m ati o n) - c/c 0 (A D E) 



Time (days) 

Fig. 3: Observed and Fitted BTC using ADE 

Fig. 3 shows the curves fitted using Microsoft Excel. 
The average pore water velocity was assumed to be a 
constant of 293.79 cm d" 1 throughout the soil profile. 
The best fit between the two curves were obtained 
when D e value was set to 500 cm d' . The 
longitudinal dispersivity (X) was calculated using Eq. 
9. The value of dispersivity (X) for bromide at the 
depth of 89 cm for the proposed sandy loam was 
calculated as 1.70 cm. This value lies well within the 
range (0.1-2 cm) reported for homogeneous saturated 
soils in laboratory column experiments (Bear, 1972). 


Summary and Conclusion 

The determination of various solute transport 
parameters is important mainly in the groundwater 
pollution studies. In present study, the dispersivity of 
bromide ions in a laboratory soil colu mn was 
determined by the inversion of Advection-Dispersion 
Equation. The dispersivity of the bromide ions in a 
homogeneous sandy loam soil colu mn was found to 
be 1.7 cm. 
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